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Ginsparg- Wilson Games 
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I implement a set of tricks for constructing lattice fermion actions which approximately realize the Ginsparg 
Wilson relation, with very promising results from simulations. 



It might be useful to have a simple lattice 
fermion action S = rpDip which approximately 
obeys the Ginsparg- Wilson (GW) relation |Q] 
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Published algorithms cost (apparently) hun- 
dreds of times as much as the usual clover action. 
I describe an approach which costs about a factor 
of 6.5 X (A*" -I- 1) as much as the clover action for an 
A''th order approximation, and even = 1 looks 
quite promising. 

The ideas in this work are based on three 
remarkable formulas first published by Neu- 
berger|^: Introducing a zeroth-order Dirac op- 
erator Do and defining z = 1 — Do/ro, a GW 
action (with R = r^) is 
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The inverse square root is approximated by 
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(3) 



(cj = cos2(7r(j + l/2)/(2A^)), = 1 - c^) and 
^ A B-A 



B B 

Here B — A — PU^^) is the polynomial 

N 

\\_(cjZ^z + Sj) - {z/N)^W{c,z^z + Si). 
To use this for propagators, note Dq^iP = < 

^ = (DW)-V = i?(w^^^))-V. 
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(i.e. "0 is found by inverting the simple differential 
operator W, and then multiplying by the local 
operator B.) Of course, one needs a Dq for which 
Eqn. H works well for small A^. 

A good -Do should already be very chiral. This 
immediately suggests that we begin with a fat 
link action-these actions are already quite chiral 
as shown by their small mass renormalization and 

Za^I) §■ 

The eigenvalues of a GW action lie on a circle. 
I determine the best Dq by taking a free field test 
action and varying its parameterization to opti- 
mize its eigenvalue spectrum (in the least-squares 
sense) for circularity, for some vq (which can also 
be varied; the optimal value is about 1.6). The ac- 
tion of choice is "planar:" it has scalar and vector 
couplings S = J2x,r i^ix)iXir) + i-ff^Pf_,{r))ip{x + r) 
for r connecting nearest neighbors (r — ±/t; 
A = Ai = —0.170, Pf_i — —0.177) and diagonal 
neighbors {r — i/t ± i/, ly ^ p; X = X2 = —0.061, 
Pf, = -0.053; A(r = 0) = -8A1 - 24A2). The 
approach of the eigenvalues to a circle is shown 
in Fig. |l]. The massive action for bare mass m 
is obtained from the m = one by|^ D{m) = 
(1 + am/2)Do + am. 

One might think that the iteration could be 
done starting with the Wilson or clover action. 
The trick of Eqn. ^ docs rapidly pull the eigen- 
modes onto a circle, but the problem is the de- 
composition into the W/B form. Unlike for the 
planar action, the eigenvalues of W are thrown far 
out into the complex plane. This is shown in Fig. 
H Unfortunately, W is the matrix which is to be 
inverted for propagators. Since the high momen- 
tum modes of a fat link action don't see the gauge 
fields very well, they behave like free field modes. 
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The wide spread of eigenvalues means that in real 
simulations, W^-^^ becomes ill-conditioned even 
for small TV. 

Chiral properties of the action, in four dimen- 
sions, are tested first by computing the value of 
its smallest real eigenvalue A on a set of isolated 
instanton configurations (the instanton radius is 
p) (Fig. H). In an exact GW action the real 
eigenvalue would be zero until the instanton fell 
through the lattice, when it would disappear. In 
an ordinary action, A is a smooth function of p, 
close to zero for big p and moving away from zero, 
generally to a positive value, as p decreases, until 
the eigenvalue collides with a doubler and goes 
imaginary. For a better action, A keeps closer to 
zero and breaks away more steeply, with a step 
function for A as the desired limiting result. 



Figure 1. Free field spectrum of d'^^, with the 
planar action as the kernel, for TV = 0, 1, 2, 3. 
Only the Im A > eigenvalues are shown. 
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Figure 2. Free field spectrum of the operator 
W'^^\ with the Wilson action as the kernel, for 
TV = 1, 2, 3. 
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Figure 3. Real eigenvalue spectrum of on 
background instanton configurations, for TV = 0, 
1, 2. 



As shown by the pion mass in Fig. ^ (quenched, 
for SU{3), a = 0.2 fm, 8^ x 24 lattice),the zeroth 
order action is already very chiral and TV = 1 
iteration is even more so. 
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Figure 4. Pion mass and quark mass (from the 
PCAC relation) for iV = 0, 1. 



In Fig. H I show the N/p mass ratio at ir/p = 
0.7 for the = and 1 versions of this action, 
along with other actions. Both of the new actions 
(on improved background gauge configurations) 
have small scaling violations for this observable. 

As an added feature, the hadron dispersion 
relation for these actions is better than for the 
clover action (the extra terms in the planar ac- 
tion can be tuned to optimize this[^). 

Finally, a rough calculation of Zj^ from Ward 
identities produces a value quite close to unity-as 
the fat link clover action |Q gave. 

To conclude: this is an approach towards the 
construction of a GW action in which all the cost 
is "up- front" in the evaluation of Dq, but the 
gain is that probably only a few terms (maybe 
just A^ = or 1) in the expansion of I?^^ are 
needed. One also only needs to invert the sim- 
ple (but messy) differential operator VF(^) (no 
inverse inside an inverse is needed). It would be 
very interesting to tackle the hard lattice prob- 
lem of i'ipip) along the lines of Ref. |7j with this 
approach. 

This work was supported by the US Depart- 
ment of Energy. 




Figure 5. N/p mass ratio at n/p = 0.7. A" = 
and A^ = 1 actions at a = 0.2 fm are fancy crosses, 
octagons-KS fermions, bursts are NP Csw thin 
link actions, squares fat link clover on Wilson 
gauge configurations, and the rest on improved 
gauge configurations. 
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